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Abstract 



C$ ' We consider a magnetic Schrodinger operator in a planar infinite strip with 

frequently and non-periodically alternating Dirichlet and Robin boundary con- 
ditions. Assuming that the homogenized boundary condition is the Dirichlet or 
the Robin one, we establish the uniform resolvent convergence in various operator 
norms and we prove the estimates for the rates of convergence. It is shown that 
these estimates can be improved by using special boundary correctors. In the 
case of periodic alternation, pure Laplacian, and the homogenized Robin bound- 
ary condition, we construct two-terms asymptotics for the first band functions, as 
well as the complete asymptotics expansion (up to an exponentially small term) 
for the bottom of the band spectrum. 



o 

1 Introduction 



In the present paper we study a magnetic Schrodinger operator in an infinite planar strip 
with frequently non-periodically alternating Dirichlet and Robin boundary conditions, 
cf. fig. 1. This model was formulated first in [4] for the pure Laplacian with periodically 
alternating Dirichlet and Neumann boundary conditions. The further studies were done 
in [1], [2], [3]. The motivation for such studies as well as the reviews of previous results 
were done in all the details in [1] and [4] and here we just briefly remind it. 

The motivation is threefold and comes from the waveguide theory and the homoge- 
nization theory. In the former one of the popular models is the waveguide with windows. 
Usually the waveguide is modeled by a planar strip or layer, where an elliptic operator 
is considered. The windows are modeled by a hole cut out on the boundary of the 
waveguide and coupling it with another waveguide. In the symmetric case the hole can 
be replaced by a segment on the boundary at which one switches the type of the bound- 
ary condition. Such models in the case of one or several finite windows were studied 
by various authors, see, for instance, [5], [6], [7], [8], [9], [10], [11]. In our model the 
windows are modeled by the segments on the boundary with a general Robin condition. 
Each segment is finite, while their total number is infinite. Exactly the last fact is the 
main difference of our model in comparison with those in [5], [6], [7], [8], [9], [10], [11]. 

The second reason to study our model comes from the series of papers devoted to 
the problems in bounded domains with frequent alternation of boundary conditions. 
Not trying to cite all of them, we mention just [12], [13], [14], [15], [16], [17], [18], [19], 
[20], [21], [22], [23], [24]. So, it was natural to study the case of an unbounded domain. 
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Figure 1: Waveguide with non-pcriodically alternating boundary conditions 



The third and the main part of the motivation is the recent series of papers by 
M.Sh. Birman, T.A. Suslina on one hand and by V.V. Zhikov and S.E. Pastukhova on 
the other - see, for instance, [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35] and 
further papers of these authors. In these papers the homogenization of the differential 
operators with fast oscillating coefficients in unbounded domains was studied. The 
operators were considered as unbounded ones in appropriate Hilbert spaces. The main 
result was the proof of the uniform resolvent convergence as well as the leading terms 
of the asymptotics for the perturbed resolvents in the sense of the operator norms. The 
last cited series considered the problems lying in the intersection of the homogenization 
theory and the spectral theory of unbounded operators. The approaches of both the 
theories were involved and it happened to be fruitful and interesting. From this point 
of view our model continues the ideas of papers by M.Sh. Birman, T.A. Suslina, V.V. 
Zhikov, S.E. Pastukhova, since we consider the problem from operator point view, but 
with the boundary geometric perturbation from the homogenization theory. 

In comparisons with our previous works [1] , [2] , [3] , [4] , the present one has several 
advantages. While in the cited papers we considered the pure Laplacian only, now 
we deal with a much more general operator involving variable coefficients. We also 
consider alternating Dirichlct and Robin boundary conditions instead of Dirichlct and 
Neumann ones in [1], [2], [3], [4]. Moreover, now the homogenized operator involves 
Robin boundary condition instead of Dirichlet or Neumann ones. More precisely, we 
have an additional term in the coefficient in the homogenized Robin condition, and this 
term appears due to the geometric structure of the perturbation. Such situation was not 
considered before and it is the most complicated among the possible ones. One more 
important difference is that in the present paper we consider non-periodic alternation 
of boundary conditions, while in [1], [2], [3], [4] only strictly periodic alternations were 
studied. Our assumptions for the structure of alternations arc quite weak and include 
quite a wide class of possible non-periodic alternations. We prove the uniform resolvent 
convergence in such case when the homogenized condition is the Robin one and the same 
result for the homogenized Dirichlet condition. The estimates for the rate of convergence 
are given in two possible operator norms, see Theorems 2.1, 2.2. We also find that in 
the case of the homogenized Robin condition the uniform resolvent convergence does 
not hold in the sense of the norm of the operators acting from L2 to W\ , while it was 
true in the case of homogenized Dirichlet and Neumann condition in [1], [2], [3], [4]. 

To obtain the described results, we use the combination of the techniques employed 
in [1], [4] for periodic case and that from [13]. Exactly in the last paper the non- 
periodic alternation in the bounded domain was treated, and we adapt these ideas 
for our case of the unbounded domain. In [13] the main result were the asymptotic 
estimates and the asymptotic expansions for the eigenvalues of the pure Laplacian, while 
the resolvent convergence was not considered. It should be also said that it is the first 
time in the problems on boundary homogenization that for non-periodic alternation the 
uniform resolvent convergence is proven and the estimates for the rates of convergence 
are provided. 

The next part of our main result concerns the description of the asymptotic behavior 
of the spectrum. Since in the case of general operator with non-periodic alternation 
the structure of the spectrum can be very complicated, we restrict ourselves to the 
case of pure Laplacian with periodic alternation. In this case we can apply Floquct- 



Bloch decomposition and get the band spectrum. Our main results are two-terms 
asymptotics for the first band functions and the complete asymptotic expansion (up 
to an exponentially small term) for the bottom of the spectrum. Similar results were 
obtained in [1], [2], [3], [4] for the case of homogenized Dirichlet or Neumann condition. 
In the present paper we prove them in the case of homogenized Robin condition assuming 
as above that the geometric structure of the alternation generates additional term in 
the homogenized condition. 

In conclusion of this section, let us describe briefly the structure of the paper. In 
the next section we formulate the problem and the main results. We also discuss open 
problems, and, in particular, the Bethe-Sommerfeld conjecture for our model, as well 
as possible ways of treating them. In the third section we prove the uniform resolvent 
convergence in the case of the homogenized Robin condition and establish the estimates 
for the rates of convergence. Similar results but for the homogenized Dirichlet condition 
are proven in the fourth section. In the fifth section we study the band functions in the 
periodic case and we prove two-terms asymptotics for them. In the last sixth section we 
construct the complete asymptotic expansion for the bottom of the essential spectrum 
in the case of periodic alternation and the homogenized Robin condition. 

2 Formulation of the problem and the main result 

Let x = {x\,X2) be the Cartesian coordinates in R 2 , Q, :— {x : < x 2 < n} be an 
infinite strip of the width n, e be a small positive parameter. By r + and T_ we denote 
respectively the upper and lower boundary of 57. On the lower boundary we introduce 
an infinite set of the points (s|,0), j S Z, sg = 0; s j < s j+i sucn that the distance 
between each two neighboring points is of order 0(1). By od~(e), a~(e), j £ Z, we 
denote two sets of functions such that 

*j-i+«j"-i( e )<*i-«7( e )> -?' eZ - (2- 1 ) 

Employing the introduced points and functions, we partition the lower boundary T_ as 
follows, 

7 e := U{x:£s s j -eaj(e)<x 1 <es e j +eaf(e),X2=0,j(='Z}, r e :=T-\^. 

j=—oo 

In this paper we study a magnetic Schrodinger operator in £7 subject to the Dirichlet 
boundary condition on T + U 7 e and to the Robin condition on T £ . We define the 
differential expression corresponding to the operator as 

H e := (iV + A) 2 + V in L 2 (0), (2.2) 

where A = A(x) = (A 1 (x),A 2 (x)) is a magnetic potential, and V = V(x) is an electric 
potential. We assume that Ai e W^(O), i = 1,2, V £ Loo(O) and these functions are 
real-valued. The boundary conditions for the operator l-i e are introduced as 

u = {) on r + U7 e , (-- L4 2 + &Wi = on T e , 

where A 2 is supposed to be A 2 (xi,0), and b € W^R) is a real- valued function. Rig- 
orously we introduce H e as the self-adjoint operator in L 2 {Q) associated with the sym- 
metric lower-semibounded closed scsquilinear form 

f} e [u,v] := {(iV+A)u,(iV+A)v) L<n) + (Vu,v) L2{n) +{bu,v) L2{re) on K^fi, r+U 7e ). 

(2.3) 
Hereinafter we denote by W 20 (Q, S) the set of the functions in W 2 (Q) having the zero 
trace on the manifold S lying in the closure of the domain Q. Our main aim is to study 
the asymptotic behavior of the resolvent and the spectrum of H e as e — ► +0. 



Before presenting the main results, we make certain assumptions on the structure 
of the alternation of the boundary conditions in the operator H e - The first assumption 
describes the distribution of the parts of j £ . 

(Al). There exists a function d E = # e (s) <E C 3 (R) such that 

e (e*|) = sttj, (2.4) 

q 1 < j?' E (s) s: ci, seK (2.5) 

K'(*)| + K'(«)l < Ci, seK, (2.6) 

where ci is a positive constant independent of s and e. 

The second assumption concerns the lengths of the parts of "f e . 

(A2). There exists a strictly positive function r\ = 77(e) and a positive constant c 2 Sj 1 
independent of e and j such that 

2c 2 r)(e) < 07(e) + ot(e) ^ 277(e). (2.7) 

Our first main result concerns the resolvent convergence of H e - We consider two 
main cases relating to various possible homogenized operators. In the first case we 
assume that 

— — -> -A", e -> +0, A" = const > 0. (2.8) 

em 77(e) 

We shall show that in this case the homogenized boundary condition on T_ is the Robin 
condition with certain coefficients. In order to formulate precisely this result we first 
introduce additional notations. 

Given a constant fi ^ 0, by H^ we indicate the self-adjoint operator associated 
with the symmetric lower-semibounded closed sesquilinear form 

&r } [«. v] :=((iV + A)u, (iV + A)v) L2(Q) + (Vu, v) L . 2{n) 

+ (bu,v) L2(r _ ) +((K + ^' e u,v) L2( r_ ) on Wl (n,T + ). 

In the same way as in [36, Sec. 3] one can check that the domain of H^ consists of the 
functions in W% (0) satisfying the boundary conditions 

u = on r+, (-- iA 2 + b+(K + n)'&' e \u = on r_, (2.9) 

while the action of the operator Hj^ is described by the same differential expression as 
in (2.2). 

By || • \\z!-+z 2 we indicate the norm of an operator acting from a Banach space Z\ 
to a Banach space Z 2 . 

Now we are ready to formulate our first main result. 

Theorem 2.1. Suppose (Al), (A2), and (2.8). Then the inequalities 

|| (H £ + i)- 1 - (H^ + irii^nj-n^n) < Ce(K + M )| ln(A + M )e|, (2.10) 

|| (H £ + i)- 1 - (H^ + ir'him^m ^ C(e(K + M )| ln(A + M )e| + u), (2.11) 



I (H £ + i)- 1 - (4 0) + iy'h^-rwim < C{K + y,) 1 ' 2 , (2.12) 



hold true, where the constants C are independent of e, and pL = /i(e) is given by the 
formula 

M = /^):=-^-*, ^^7^+0. (2.13) 



There exists a boundary corrector W = W(x,e, fi) defined explicitly in (3.24) such that 
the inequality 

\\(H £ + i)- 1 - (1 + W)(«&° + ir'h^a^wiia) < Ce(K + /x)| \n(K + »)e\, (2.14) 

holds true, where the constant C is independent of e. 

Let us discuss the results of this theorem. First of all we note that the operator 
"Hpf* depends on e not only via /i, but also due to the presence of the boundary term 
(K + fi)-ff' e in (2.9). However, wc can easily get rid of such dependence, if we additionally 
assume that 

there exists a bounded function $» : R — >• E, such that sup \&' e — $*| — > 0. (2.15) 

R 

Then we can replace the mentioned term by {K + (1)$*, and it will also generate addi- 
tional term C(K+n) sup \&' e — $*| in the right hand sides of the estimates in Theorem 2.1. 

We track the dependence of K in the right hand sides in the estimates in Theorem 2.1 
just to involve also the case K = 0, where, as we see, the rates of the convergence 
are better. We observe that the estimate (2.12) is of interest only for K = 0, while 
otherwise it states no convergence. Moreover, as K ^ 0, there is no even strong resolvent 
convergence of % e to T-C^ in the sense of the norm || • || L 2 (n)^w 1 (si) ■ Indeed, assume for 
simplicity that $ e is independent of e and for a given / e L 2 (Q) let u e := (H e + i) -1 /, 
uq := (rift + i) -1 /. Then by the definition we have 



(iV + A)u £ \\ 2 L , n) + (Vu £ , Oi 2 (fi) + (bu e ,u e ) L2{re) -i|MI| 2(0 ) = (f, u e)L 2 (n), 



\L 2 (n) T V u>s,Us)L 2 (n) T \ua e ,a e)L2 (Y c ) - H\^s\\L 2 {0.) 

|| (iV + A)u \\l 2m + (V«o, u )l 2 (q) + {(b + K§')u Q , u ) L2{r) 

2 



i|l u o|lL 2 (n) = (/,wo)i a (n)- 



(2.16) 



If we suppose the strong resolvent convergence at least for a sequence e — > +0, we can 
pass to the limit in the former identity and we get 

j|(iV + ^)u |lL(o) + ( Vu 0' u n)L 2 (n) + (bu , uo)i 2 (r E ) -i|N||£ 2 (n) = (i>o)L 2 (n)- 
This contradicts to (2.16). Nevertheless, due to (2.11) the uniform resolvent convergence 
holds in the sense of the norm || • ||i 9 (n)->L9(0) no matter whether K vanishes or not. 

We can see also from Theorem 2.1 that the estimates for the rates of convergence 
depend highly on the operator norm and on the approximating operators. The worst 
estimate is in (2.12). The first way to improve the convergence is to replace the norm 
by a weaker one as it was done in (2.10) and (2.11). The former estimate is better than 
the latter, but the price to pay is the dependence of the approximating operator rt\^ 
on fi. At the same time, this dependence is quite simple and it is only in the boundary 
condition (2.9). If we keep the operator norm || • \\L 2 (n)->w 1 <n)i we have to employ a 
special boundary corrector W in order to get the reasonable estimate, see (2.14). Here 
the rate of the convergence is the best one among given, but the approximating operator 
is the most complicated. We also note that these effects for pure Laplacian and strictly 
periodic alternation were already found in [1] in the particular case K = b = 0. 

In the second case we have K = 00, or, more precisely, 

eln?7(e)^-0, e -> +0, (2.17) 

and then the homogenized boundary condition on T_ is the Dirichlct one. In this case 
the homogenized operator is denoted as %d and is introduced as the self-adjoint one 
associated with the symmetric lower-semiboundcd closed scsquilinear form 

^D[u,t;]:=((iV + A)«,(iV + i4)«) ia(n) + (Vu,t;)i a (n) on ^ (fi,r_). 

Again by analogy with [36, Sec. 3] one can check that the domain of %d consists of 
the functions in W$ (O) vanishing on <9£7, while the action is given by the differential 
expression in (2.9). 

In this case we replace the assumption (A2) by a stronger one. 



(A3). There exists a strictly positive function 77 = 77(e) and a positive constant C3 inde- 
pendent of e and j such that 

77(e) <af(e) < c 3 n(e) s^ -. 

The result on the uniform resolvent convergence is formulated in 
Theorem 2.2. Suppose (Al), (A3), and (2.17). Then the inequality 

UHt+i)- 1 - (Ho +ir 1 ||L 2 (fi)^(n) < Ce^dlnsin^e)! +cos77(e)) 1/4 , 
holds true, where the constant C is independent of e 

With respect to Theorem 2.1, in Theorem 2.2 the estimate for the rate of conver- 
gence is given in the best possible norm and no corrector is involved. Here the rate of 
convergence can be improved only by assuming 77(e) — > ir/2 — 0, since in this case the 
function | In sin 77(e)! + cos 77(e) tends to zero. This is quite natural, since the assump- 
tion 77 = 7r/2 means r e — and r y £ = r_. Then there is no alternation of boundary 
conditions in this case and H e = Hv ■ 

The proven theorems and [37, Ch. VIII, Sec. 7, Ths. VIII.23, VIII.24] imply the 
convergence for the spectrum and the spectral projectors ofH e . By er(-) we denote the 
spectrum of an operator. 

Theorem 2.3. The spectrum ofH e converges to that of Ho, where 

H Q := H { £\if we assume (Al), (A2), (2.8), (2.15), and (2.9) holds for #* in place of &' e , 

or 

Ho '■= Hn, if we assume (Al), (A3), and (2.17). 

Namely, if X ^ a(Ho), then X ^ o~(H e ) for e small enough. If X £ a (Ho), then there 
exists X 6 £ o~(H £ ) such that X e —)■ A as e — ► +0. The convergence of the spectral 
projectors associated with H e and Ho 

WP(a,b)(H e ) - ■p(a,b)(H )\\L 2 (n)^L 2 (Q) -+0, £ -> 0, 

is valid for a,b g" o~(Ho)> a < b. 

The next part of our results concerns the case of the pure Laplacian with the periodic 
alternation, namely, 

A = 0, V = 0, 6 = const, H e = -A, t? e (s) = s, af(e) = r)(e). (2.18) 

In this case the operator H e is periodic due to the periodicity of the boundary conditions. 
We apply the Floquet-Bloch expansion to characterize its spectrum. We consider the 
cell of periodicity f2 e := {x : \xi\ < stt/2, < x 2 < it}, and let -% := dft e n j e , 
r e := dVt e n r e , T± := <9fL n T±. By H e (r) we denote the self-adjoint operator in 
L 2 (£l E ) associated with the lowcr-semiboundcd closed symmetric sesquilinear form 

Ur)[u,v] := ((iA - I) u, {l±- - I) «)^+(|£. §xk) L2{n / bM ^ 

on W\ per (fl £ , r + U7 e ), where r £ [— 1, 1) is the (rescaled) quasimomentum. The symbol 
W 2 1 , P er-(^ejr + U %) indicates the set of the functions in W 2 1 (f2, r ,r+ U %) satisfying 
periodic boundary conditions on the lateral boundaries of f2 e . The operator H e (t) has a 
compact resolvent and its spectrum consists of countably many discrete eigenvalues. We 
denote them as X n (r, e) and arrange in the ascending order counting the multiplicities. 
By [4, Lm. 4.1] we know that 

oo 

a(H E ) = a c (H e ) = \J {A„(r, e) : r £ [-1, 1)}, 



where er c (-) indicates the essential spectrum of an operator. 

Let £ £ be the subspace of L 2 (^l e ) consisting of the functions independent of x\. We 
decompose the space L2(£l e ) as 

L 2 (n s )=£ s ®£^, 

where £^r is the orthogonal complement to £, E in L 2 (£l £ ). By Q M we denote the 
self-adjoint operator in £ e associated with the lower-semiboundcd symmetric closed 
sesquilinear form 

qii [u,v]:=(-^,-p-) + (b + K + f x)u{0)^) on W^ ((0,7r),M). 

\dx 2 dx 2 J L2i0 ^ ) 

An alternative definition of Q^ is as the operator —-j^i in L2(0,tt), whose domain 
consists of the functions in W 2 (0,ir) satisfying the boundary conditions 

u(tt)=0, u'(0)-(b + K + n)u{0) = 0. 

The uniform resolvent convergence for "H £ (r) and the asymptotic behavior of the 
band functions A„(r, e) read as follows. 

Theorem 2.4. Suppose (Al), (A2), (2.8), and (2.18). Let \t\ < 1-x, where < x < 1 
is a fixed constant. Then for sufficiently small e the estimate 

T 2\-i 
n e (T)-—\ -Q~ 1 ®0 ^Ce 1 ' 2 {K + ^){x- 1 / 2 ^\\ne{K + ^)\), 

(2.19) 
holds true, where /i = fi(e) is defined by (2.13), and C is a constant independent of 
e. Given any N, for e < 2m: 1 ' 2 N^ 1 the eigenvalues A„(r, e), n = 1, . . . , N, satisfy the 
asymptotics 



T 2 



A„(r, e) = — + A n (fi) + R n (r, e, /x), 2Q 

\Rn{r,e, n)\ ^ Cn A s x ' 2 {K + M )(^ 1/2 + | lne(A' + fj)\), 

where A„(/i), n = 1, ... , N, are the first N eigenvalues of Q^, and the constant C is 
the same as in (2.19). The eigenvalues A n (/x) are holomorphic w.r.t. /i and solve the 
equation 

a/Acosn/A7T+ (K + b + /J,) sin VXtt = 0. (2.21) 

The proven theorem implies immediately 

Corollary 1. Suppose (Al), (A2), (2.8), and (2.18). Then the length of the first band 
in the essential spectrum ofH e is at least of order 0(e~ 2 ). 

We mention that similar results for the homogenized Dirichlct and Neumann condi- 
tion were established in [4, Sec. 2], [1, Sec. 2]. 

In addition to the asymptotics of the first band functions, we describe the asymp- 
totics for the bottom of the essential spectrum of H e . We let 

0(ti,ta):=2ti53 



' nUAri 2 ~t 2 + 2n + t 1 ) 

oo (2-22) 



£-[ n(^An 2 -t 2 + t 1 ) (V^ri 2 - t 2 + 2n + t x ) 
We shall prove in Lemma 6.2 that the function 9 is jointly holomorphic in t\ and £2- 



Theorem 2.5. Suppose (Al), (A2), (2.8), and (2.18). Then the identity 

inf Ai(r,e) = Ai(0,e) (2.23) 

re[-l,l) 

holds true. The bottom Ai(0,e) of the essential spectrum a c {T-L e ) has the asymptotic 
expansion 

A 1 (0,£) = A( £ , M (e))+o((/v+ M )c- 2£ "+r / 1 / 2 (X + / i)+ £ 1 /V /2 (^ + M) 1/2 ), (2-24) 
where A(e,/x) is the root of the equation 



A cos VAtt + (b + K + /i) sin y/Xn) (l - e(K + ^)6(eb, e 2 A)) 

+ e(K + fi) 2 6(eb, e 2 A) sin VAtt = 0, 



(2.25) 



satisfying 

A(e,fj l )=A 1 (e,fj l ) + o(l), e -> +0, ^ -> +0. (2.26) 

The function A is jointly holomorphic in e and /i and can be represented as the uniformly 
convergent series 

oo 

A(e,„) = A 1 (i J ,)+J2^M, (2-27) 

where the functions Tj are holomorphic in [i. All these functions can be found explicitly 
and, in particular, we have 



TiG«) = 

t 2 (m) = 



7T 



Ai(a*)(Ai(a*)+m)' 



6 n(K + 6 + ^) 2 + Ai7r + 6 + K + (i ' 
C(3) A 1 (/.)(A' + M ) 2 (A 1 (,i)+2& 2 ) 



(2.28) 



4 tt(A' + b + /i) 2 + Aitt + b + K + ju' 
where £(t) is the Riemann zeta- function. 

The first result of this theorem, namely the identity (2.23), was proven in [1, Sec. 2] 
and [4, Sec. 2] for the case of the homogenized Dirichlet and Neumann condition. The 
main result of Theorem 2.5 is the asymptotics (2.24). Although it contains just one 
term A{e,jj), this term is holomorphic w.r.t. e and fi and is defined quite explicitly 
as the solution to the equation (2.25), (2.26). This solution can be represented as 
the convergent series (2.27) and the latter can be regarded as the asymptotics for A. 
In other words, the power in s and 77 part in the asymptotics for A can be summed 
up to the exponentially small error, see (2.24). It is a strong result for the problem 
in the homogenization theory, since usually complete asymptotic expansions in the 
homogenization theory can not be summed up. We also mention that previously similar 
result was only known in the case K = b = 0, sec [1, Sec. 2]. 

In conclusion to this section let us discuss some open questions related to the studied 
model. The main question is the gap conjecture, namely, the question on the gaps in 
the essential spectrum of r H. e . It is not known whether such gaps exist or not. If exist, it 
would be interesting to prove the Bethe-Sommcrfcld conjecture for our model, namely, 
that the number of the gaps is finite. The next question would be the dependence of the 
number, the location, and the sizes of the gaps on e. These questions were formulated 
even in the first paper [4], where the studied model was proposed. Our conjecture was 
that the gaps exist and a possible way of proving it was to study in more details the 
behavior of the band functions as e — > +0. In fact, the idea was to extend the theorems 
like Theorem 2.4, trying to include into consideration vicinities of the points t = ± — 1, 
to construct the complete asymptotic expansions for A„, and to get these results not 
only for the first ones but for all band functions A„. Now we have to admit that these 
questions are much more complicated as we expected in the beginning. In particular, 



one needs to develop a new technique in comparison with ours to construct the complete 
asymptotic expansions and/or to include the vicinities of the points r = ±1. Moreover, 
there is a usual problem with ordering of the band functions. Indeed, even in the case 
b = K = 0, if we consider the homogenized operator as the periodic one with the period 
en, then the associated band functions are ^ m + T > — 1_ ( n _|_ i^ ; TO) n g ;g + _ Then the 
issue of putting these eigenvalues in the ascending order for all values of e is not trivial. 
Attempts to improve (2.19) seem to be not useful, since such convergence can not 
describe in an appropriate way the desired asymptotic behavior of the band functions. 
This is why to solve the gap conjecture for our model, one has to develop a new approach. 
It is likely that the approaches used in the proofs of the Bethe-Sommerfcld conjecture 
for various operators can be useful. Moreover, now we conjecture that our model has no 
gaps in the essential spectrum. The first argument supporting such conjecture is that 
the homogenized operator has no gaps at all and the number of overlapping bands in 
its essential spectrum increases as the point goes to infinity. The second argument is 
that we can make the rescaling in the original operator x *-¥ xe . Under such rescaling 
the operator Tt e becomes the operator —A in the strip {x : < x 2 < ne^ 1 } with 
fixed alternation of the boundary conditions. As e — > +0, the new strip "tends" to the 
half-space x 2 > 0. Hence, one can expect that the spectrum of the rescaled operator 
converges in certain sense to that of — A in x 2 > with the original fixed alternation. 
And since the spectrum of the latter is always [0, +00), we can conjecture that there is 
no gaps in the perturbed spectrum. If the essential spectrum indeed has no gaps, then 
Theorem 2.5 and its analog in [4] gives the complete description of the location of the 
essential spectrum, namely, of its bottom. 

3 Resolvent convergence: homogenized Robin condi- 
tion 

In this section we prove Theorem 2.1. 
Given a function / € £2(0), we denote 

u e := (H e - i)" 1 /, u^ := (H^ - i)" 1 /- 

We first prove an auxiliary lemma which will be one of the key ingredients in the 
proof of Theorem 2.1. 

Lemma 3.1. Let W = W(x, £,/x) be a real function belonging to 

C*(Ii)nC 2 (n\{x:i 1 =es e j ±eaf(e), x 2 = 0, j 6 Z}), (3.1) 

satisfying boundary conditions 

dW 

W = -\ on le , =-(K + ii){l + e Ve )$' E on T e , (3.2) 

0x2 

bounded uniformly in CI, and having differentiate asymptotics 

W (x, s, fi) = cf{e, n) yfff sin -L + 0(rf), rf -► +0. (3.3) 

Here ip e = ip 6 (xi) € C(R) is a bounded function, (r- , #• ) are polar coordinates centered 
at (es^ ± ea (e), 0) such that the values $■ = correspond to the points of ~f e . Assume 
also that AW € C(fi) and the function AW is bounded uniformly in fi. Let 

v £ :=u e -(l + W)u { >*\ (3.4) 



Then (1 + W)u^ belongs to W^ )0 (O,r + U j e ), and 

||(V + L4)t> £ ||! 2(n) + (Vv e ,v e ) L2m +(bv E ,v e ) L2( r e ) +i||w e ||i 3 (n) 

=(f,VeW) L2m + (u^AW,v 6 ) L2in) ~ 2i(u^W,v e ) L2in) 

- 2{Vu^\Wv E ) L2(n) - 2(W(V + iA)u^\ (V + iA)v 6 ) L2{n) ( 3 ' 5 ) 

- ((2b+(K + iJ,)(l + eip £ )-d' e )u M ,Wv £ ) 

Proof. The functions u e and w^ satisfy the integral identities 
((V + L4)u £ ,(V + iA)<t>) L2m + {Vu 6 , <j>) L2iQ) 

+ (bu E ,ct)) L2{re) +i(u e ,0) L2 (O) = (/,0)l 2 (O) 
for all $ e Wlo(Cl, r+ U 7e ), and 

((V + L4)«<"\(V + iA)4>) L2(n) + {Vu^^) L2{n) 



(3.6) 



+ ((b + (K + l x)d' e )u^\ct>). ) +i(u^\4>) L2m = (/,0) i2( n) 



(3.7) 



for all </> € Wj o(^> r+). By analogy with the proof of Lemma 3.2 in [4], one can show 
that (1 + W)4> G W2 (f2, T + U 7 e ) for each <p which belongs to the domain of either T-L £ 

or w£' . Hence 

(1 + W)u w G W 2 1 ,o( n > r + U 7£ ), (1 + W)v e G W r a 1 I o(fi.r+ U 7*)- 
We choose the test function in (3.7) as <f> = (1 + W)u e . It yields 
((V + iA)vP>, (V + L4)(l + W)v e ) L2{n) + {Vu^\ (1 + W)v £ ) L2m 

+ ((b + (K + ^' e )u^\ (1 + W)v E ) L2{r) + i(«M (1 + W0« C )£,(O) 
=(/,(l + WK)L a (0) S 

and we rewrite this identity as 

((V + iA)uP>, (1 + W)(V + i4K) La(0 ) + (Vu^\ (1 + W> e ) L2(n) 
+ (bu^\ (1 + W> e ) L2(r _) + i(u (M) , (1 + W> e ) L2(n ) 
=(/, (1 + W)v 6 ) L2{n) - ((V + iA)u^\v e VW) L2{a) 

and 

((V + L4)(l + W)vP\ (V + L4H) i2(n) + (V(l + W> (m) ,Omh) 

+ (6(1 + W)u w ,t; e )L 2 (r_) + i((l + W)u^\v e ) Laia) 
=(/, (1 + W> e ) La(0 ) - ((A" + m)^ M , (1 + ^)« £ ) L2(re) 
- ((V + iA)u^\v e VW) L2{n) + (« M VWi (V + L4K) La( n). 
We deduct the last identity from (3.6) with <j> = v 6 and get, 

||(V + L4H||! 2(n) + {Vv e ,v e ) L2{n ) + {bv e ,v e ) L2 (T s ) + i|KHi a (n) 

= - (/, W« e ) ia( n) + ((if + ^ e uM, (1 + W)v e ) L2{Te) (3.8) 

+ ((V + \A)u^\v e VW) L2(n) - (yP>VW, (V + iA)v £ ) L2(n) . 

Now we consider separately the two last terms in the right hand side of equation (3.9). 
We integrate by parts and employ the boundary condition (3.2), the identity (3.7), and 
the belongings (3.4), 

{(V + \A)u^\v e VW) L2(n) ~ (« (fl) VW,(V+iiK) i!(sl) 
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= ((V + iA)uM,v E VW) La{ n) + [uM^Vedx! 

+ (diYU M VW,v e ) L2ia) - {u^VW:iAv e ) L2{n) 
= 2((V + iA)u^\v e VW) L2{n) + (<uV>AW, v 6 ) L2{n) 

and 

((V + iA)««, v e VW) L2{n) =((V + L4)«<">, (V + \A)v e W) ia(n) 

-((V + ii)«W,^(V + iiK) Mfi) 

=(/, W« e )£ a (n) - (^ W ,^)i j( o) - h> (M) , W^ e ) L2(0) 
~((&+(A' + /i)<)>),^ e ) L2(re) 

-((V + L4)«^,W(V + L4)t; e ) ia(n) . 

These equations and (3.8) imply (3.5). □ 

Our next step is to construct the function W satisfying the hypothesis of the proven 
lemma. Then we shall be able to estimate the right hand side of (3.5) and as a result 
the W% (O)-norm of v e . Exactly the last estimate will allow us to prove Theorem 2.1. 

To construct the function W, we employ the technique from [13], [14], [23], [24] based 
on the method of matching asymptotic expansions [38] , the boundary layer method [39] , 
and the multiscale method [40] . For the periodic case similar approach has already been 
successfully applied in [1, Sec. 3]. In order to treat the nonperiodic case we employ some 
ideas in [1] in combination with the technique of [13]. 

We seek W as a formal asymptotic solution to the equation 

AW = in fl 

satisfying the boundary conditions (3.2). This solution is constructed as a sum of an 
outer expansion, a boundary layer, and an inner expansion. We begin with the boundary 
layer. 

The main idea of introducing the boundary layer is to satisfy the Neumann condition 
on r e . We first define the rescaled variables 

e=(£i,6), ei = ^^, 6 = <(*!)-• (3.9) 

£ £ 

This change of variables maps the points x = (es|,0) into a periodic set £ = (ttj,0), 

while the Laplacian becomes e -2 ($' e (xi)) A^ + C(e _1 ). These two facts are exactly 
the motivation of defining the rescaled variables for the boundary layer by (3.9). We 
construct the latter as 

W?(Z,xi,ii) = e(K + ri(l + e<p e (x 1 ))X(t), (3.10) 

and the function (p s will be determined below. 

We substitute (3.10) into the boundary value problem for W, replace then the set 
T e by T_ \ |J (espO), and equate the coefficients at the leading power of e. It implies 

the boundary value problem for X, 

BX + °° 

A ? X = 0, 6>0, ar = -l> £er°:=U:6=0}\ (J (errj.O). (3.11) 



d£ 



3- 



We are interested in a solution to this problem which is £7r-periodic in x\ and decays 
exponentially as £2 — > +00. Such solution was found explicitly in [21], 

X(0 := Relnsin(6 + i£ 2 ) + In 2 - £2, (3.12) 
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where the branch of the logarithm is fixed by the requirement In 1 = 0. Some additional 
properties of this function were established. Namely, it was shown that the function X 
belongs to C°°({£ : £2 > 0} U T ) and satisfies the differentiable asymptotics 

X(0 = ln|e-(^,0)|+ln2-e 2 +O(|e-(^,0)| 2 ), £->(7rj,0), j € Z. (3.13) 

The constructed boundary layer does not satisfy the required boundary condition 
(3.2) on 7 6 and it has also logarithmic singularities at the points (es^,0) because of 
(3.13). This is the reason why in a vicinity of these points we introduce the inner 
expansion and construct it by using the method of matching asymptotic expansions. In 
a vicinity of each point (ss £ a, 0) we introduce one more rescaled variables, 



= (4 j) ,4 j) ), 4 J) ■= (£1 - «J)n-\ 4 j) ■= 6T 1 - 



-00 



It follows from (3.13), (3.10) that 

Wf^, Xl ,ti) = - 1 - ei Ps (x 1 ) + s(K + l x)(l + S i Pe (x 1 ))(ln\^\+ln2- r^ } ) 

+ 0(s(K + n)r] 2 \<; (j) \ 2 ), e-^W.O). 

(3.14) 
Hence, in accordance with the method of matching the asymptotic expansions we should 
construct the inner expansion for W as 

W^' j {6 j \xi,n) = -l + eY^(^\e), (3.15) 

and the functions Y^'{<;^\e) should satisfy the asymptotics 

y(%a) )£ ) = (K + ij,)(l+eip £ (x 1 ))(]n\^\+]n2)- i p e (x 1 ) + o(l), \^\ -► 00. (3.16) 

We substitute (3.15) into the boundary value problem for W and equate the coefficients 
at the leading orders of the small parameters e and r\. It leads us to the boundary value 
problem for Y^\ 

A fO) F« = 0, ? 2 (i) >0, 

F W) = 0, ?<*> e 7 lj := W J) : ~2aJ(e) < <;[ j) < 2a+(e), <£> = 0}, 

^ = 0, ? «er^:=o ?1 \^ (3-17) 

± J) e {es^±eaf{e))~Ue^) 
^ (£):=± fej ■ 

We observe that due to (2.1), (2.4), (2.5), (2.7) the functions a~~(e) satisfy the estimate 

\af(e)\^C, 

where the constant C is independent of e and j. 

In [21] a special solution to the problem (3.17) was found for the case ocj{e) = ±1 
being 

Y(q) :=Rcln(z + v / ^ 2 -l), z = ^+k 2 . (3.18) 

Here the branch of the square root is fixed by the requirements y/l = 1. It was shown 
that the function Y belongs to C°°({<; : <T2 ^ 0, ^7^ (±1,0)}) and satisfies the asymp- 
totics 

Y( ? )=ln| ? |+ln2 + (9(| ? |- 2 ), ?^oc. (3.19) 

By a trivial change of variables we can obtain then the solution for our case, 
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Hence, the function Y^ has the same smoothness as Y and 

y(%(j) je ) ={K + /i)(l + ecp e {xi))(ln\^\ + In 2 - Ind,-(e) 

m x ( 3 - 2 °) 

+ (07(e) - at( £ ))^L_ + 0(| ? W|- 2 )J, | ? 0')| _, +0O . 

We compare this asymptotics with (3.16) and see that the function ip e should satisfy 
the identity 

(K + fj,) (l + eip e (xi)) In d j (s)=(p e (x 1 ) as x\ - es^ = 0(erj). (3.21) 

We define such function as follows. Let \ be an infinitely diffcrentiable cut-off function 
with values in [0, 1], being one as t < 1/4 and vanishing as t > 3/4. Denote 

M £se i + £ 4( £ )) - M £s j ~ £Q 7( £ )) 

aA£) '- 2e V (e) 

and let 

g e (6):=d j+ i(e)-x(jA(d j+1 (e)-dj(e)), enj ^ 8 ^ eir(j + 1), jgZ. (3.22) 

It is easy to see that in a e-ncighborhood of each point es 6 j we have <7 E ($ E (a;i)) = dj(e). 
Then we replace dj in (3.21) by g e ($ e (xi)) and determine ip e , 

^ := i-^ + ^m^^))- (3 - 23) 

It follows from the assumptions (Al), (A2) that the functions g e ('d e (xi)) and i? c (xi) are 
bounded uniformly in xi and e. 

As we see from (3.14), (3.20), the term —e(K + /j,)^ is not matched with any 
term in the inner expansion. It was found in [14], [23], [24] that this term should be 
cither matched or canceled in order to have a good approximation by the constructed 
asymptotics. This is also the case for our problem and to compensate it, we add the 
term e(K + /j,)^ in the boundary layer and we introduce the term — e(K + /J,)x2 as the 
outer expansion. With this trick the final form for W is 

W(x, e,n)=- {K + n)(l + e(p e (x 1 ))^' s (x 1 )x 2 

+ s(K + ^(l+s^(x 1 ))(X(0+^) (1-X(|?°'V /4 )) 



(3.24) 



+OO 



+ E x(k°V /4 )(-i + ^ ( %°V))- 



This function is well-defined since the product and the sum in its definition are always 
finite. 

Let us check that the function W satisfies the assumption of Lemma 3.1. By direct 
calculations we check that the belonging (3.1), the boundary conditions (3.2), and the 
asymptotics (3.3) hold. It follows from (3.1), (3.2) that AW € C(Cl). It is also easy to 
check that this function is bounded uniformly in Q. Therefore, we can apply Lemma 3.1. 

Our next step is to estimate the right hand side of (3.5). In order to do it, we need 
several auxiliary lemmas. The first two of them were proven in [1] for the periodic case 
# e (s) = s, af(s) = 1. For our non-periodic case the proof is the same with several 
obvious minor changes. This is why we provide these lemmas without proofs. 

Denote 

+00 +00 

n s -.= (J {x:\z-(8 c jt o)\<8}nn, j 5 e -.= \J {x-. |a-^| <S, x 2 = 0}. 
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Lemma 3.2. For any S g (0,7r/2), u e W^i^l), v € W|(0) the inequalities 
\\u\\ L2m ^ C8(\]n6\^ + l)\\u\\ wm , 

hold true, where the constants C are independent of e, 8, u, and v. 

Lemma 3.3. For 5 e [3?7 1 / 4 /4, 7r/2) the estimates 

\W\^Ce(K + /j,)(\]n8\ + l), x€il\n s , \W\^C, xen s , (3.25) 

are valid, where the constants C are independent of e, [i, r\, 6, and x. 

Due to the non-periodic structure of the alternation of boundary conditions, we need 
to prove the next two lemmas. Denote 

fl sd := jx : |£i -irj\ < -J. 

Lemma 3.4. The estimate 

\\AW\\ L2{nsj) t:Ce(K + n) 
holds true, where the constant C is independent of e and j . 

Proof. Throughout the proof by C we denote inessential constants independent of e and 
j. We first calculate AW employing the equations in (3.11), (3.17) for X and Y^\ 

AW = -{K + mM(1 + e<p e )# e )" 

+ £ {K + n)(2 v 'X^-+e^{X + i 2 )\ jf (l - *(k°V /4 )) 

^ ' J — — OC 

+ 2 J2 v x (k°V /4 ) • vw-7 at + ]T V^ mat A x (k (j) |77 3/4 ), 

j— — oo j — — oo 

W™\x,e,ll) = -l + sY^\^\e) - s(K + M )(l + e <p e (xi)){X@ + £ 2 ), 
where the arguments of Y are ( (q + ctj(e) — af(s))dj 1 (e), q!f dj 1 (e) 1 . 

3 1/4 

For x £ Qi v the most part of the terms in AW vanishes and we have 

AW = e(K + M ) (V e 0*|? + W"A . x e n s . 

It follows from (2.5), (2.6) that 

\<p' e \^Ce-\ l^'KCe" 2 , (3.26) 

where the constant C is independent of e and x\. Moreover, since by (3.22), (3.23) the 
function <7 e ($ e ) is constant as |£i — nj\ < 7r/4, the derivatives of (p e vanish for such 
values of x\. Thus, for some small fixed S £ (0, 7r/2), 

M^M^no*^) * c c 2 ( K + rf\\x\\w m : i2 >o, 6<M<.m) (327) 

where the constant C is independent of e and j . 

As x <G El j \ fli 71 , the functions xOC^I 7 ? 3 ^ 4 ) are identically one and thus 



AW = -{K + n)x 2 ((l + etp e )0' e )", x e n ed \ Hi" 
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Together with (3.26), (2.6) it implies 

'M-!>:k«»i<^ 1/4 }) 



A W\\ Lo( r, cU)l<iril/il] < C{K + ,i)7// 2 . (3.28) 



3 1/4 1 1/4 

It remains to obtain the estimate in the intermediate region $7i T ' \Qi v . Here by the 
matching performed above and, in particular, by (3.10), (3.14), (3.21), (3.23), (3.24), 
the estimates 

W7 at (x, £ ) = o( £ (A' + M )(l + £ ^(^i))(ie| 2 + kr 1 )), 
Wf at (x, e) = 0((K + M )(|el + rtl~ 2 )) = 0((K + m)^) 

are valid. It implies 

||2v x (k (J) |7? 3/4 ) ■ vwr^wr^xCk^l^H^^i^-^ni,!-") < c(^ + m)V /2 , 

while the other terms in AW can be estimated in the same way as in (3.27), (3.28). 
The last estimate, (3.27), and (3.28) imply the desired estimate. □ 

Lemma 3.5. Each function u £ W% (fi) belongs to C(f2) and satisfies the estimate 



ENIcoL^^NIi 






where the constant C is independent of e and u. 

Proof. By the standard smoothness improving theorems (see, for instance, [41, Ch. I, 
Sec. 6, Th. 3]) each function u £ Wf (fi) belongs to C(f2). Moreover, we let 



Qj 



[ ($1,3:2) : \(i-nj\ < -, <x 2 < -j, "(£1,2:2) := w(aJie S^a). 



Again by [41, Ch. I, Sec. 6, Th. 3] there exists a constant C independent of e, j, and 
u such that 



lull 2 - 



< C||«||wj(<j) < C , (lNlw, a (n i _i) + INIw?(n,) + Nlw*(n i+1 : 
where Q is a domain with infinitely differentiablc boundary satisfying 



SljCQc Uj-i U Uj U fij+i. 
Since 



we get 



Mlwf(n 3 ) ^ e ll M llw|(n eiJ )' 



l u llc'(n.,o ^ Ce '(ll^ll^K^-!) + II^H^Kn^) + Nl^(n e , 3 - +1 ))- 



We sum this estimate over j 6 Z and it completes the proof. D 
Employing Lemmas 3.2, 3.3 and proceeding as in [1, Eqs. (3.26)-(3.32)], we obtain 

h (jx) \\wS(.a)<0\\f\\ La(a) , (3.29) 

|(/, Wv s ) L2{n) \ < Ce(K + n)\]ne(K + p)|||/|U a( n}IHIw*(n), (3-30) 

\\u M W\\ L2{n) + ||WVu^|| La( n) ^ Ce(K + M )| lne(A + M)|||/||L a (n), (3-31) 

I («(%, « £ )l 2( o)+(W(V + L4)u<"\ (V + L4)« 8 ) ia(n) I (3 32) 
sC Ce(A + /u)| lne(A + /i)|||/|U 2 (n)||%||wi(n)- 
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The inequalities (3.29), (3.30), (3.31) also imply 

\\Wv e \\ L2{n) sc Cs(K + n)\]ns(K + u)\\\f\\ L2m \\v e \\ wi(Q) , (3.33) 

\(Vu^\Wv e ) L2{n) \ < Ce(K + fj,)\]ns(K + u,)\\\f\\ L2m \\v e \\ L2m . (3.34) 

To estimate the term ||u^^AM^||i 2 (n) we employ Lemmas 3.4, 3.5 and (3.29), 

\\uMAW\\ L2m = lj2h MA W\\l 2(ae A 

(3.35) 

< Ce^\K + v)\\u^\\ wm ^ Ce x l\K + M)ll/IU a (n). 

We employ Lemmas 3.2, 3.3 with 8 = e 2 (K + fi) 2 , the embedding of W£{Q) into L 2 {T-) 
and obtain 




e^CK + ^H 



(3.36) 



^Ce\K + ^\\ne{K + ^\\u^\\l 2{Te) 



+ Ce 2 (K + rf\\u^\\ 2 wim 
^Ce 2 (K + rf\lne(K + ri\ 2 \\f\\ 2 L2{n} . 

Hence, by (Al), (3.23), and the boundedness of b, we have 

((2b+(K + ri(l + Et p' e )# e )uM,Wv e ) T /r x < h (M) W|U 2( r e )|ke|| i2( r_) 
< C||^V|U 2( r e) 11^11^(0) < Ce(if + n)\lne(K + M)lll/IU a (n)IMw2(0)- 

We substitute the last estimate and (3.30), (3.31), (3.32), (3.33), (3.34), (3.35) into 
(3.5). Then we take the real and imaginary parts of the obtained estimate and we use 
the obvious inequality 

l|V«||! 2( o) < C(||(V + L4)«||i a(0) + |Mli a(n) ) (3.37) 

which is valid for any v £ W^i^l). As the result we get 

\\ v e\\l 2 (n) < Ce ( K + ^)l ln£ ( A " + M)l!l/llL 2 (0)ll«e||w 2 i(0), 

||V«e||i a( n) < Ce{K + n)\hxe(K + A0lll/IU a (O)Kllw2(n) + C||Vt; e ||i a(n) . 

Thus, 

WveWwm < Ce(tf + M )| lne(X + M)lll/IU 2 (n), 

and it proves (2.14). 

In order to prove the other estimates we shall make use of one more auxiliary lemma. 

Lemma 3.6. The estimate 

l|V(u<">W)|U 3( n) < C(K + m) 1/2 ||/IIl 2(S2 ) 
holds true. 
Proof. We integrate by parts taking into consideration (2.9), and (3.2) 

\\W(u^W)\\ 2 r Im = - ( ^-u M W,u M w] - (a(u m W),u m w) 

2{ ' \OX 2 /i 2 (r_) V /L 2 (Q) 
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L 2 (-te) 



u w!?w m w 



\ dx 2 /Lo(T_-) V ' 



Mr-) 
u M AW,u M W] -2(wVu {ti \u M Vw) 

/L 2 (Q) V J L 2 



M^) 

L 2 (fi) 
(O) 



-2(T^Vm (/j) ,u ( ' j) VW) 
£ 2 (fi) V /l 2 (o) 



- ((-L4 2 + 6 + K + //)u (,l) W, w (M) ^) L9(r 

- (u M AW,u M W 
We take the real part of this identity 

\\V(u M W)\\l 2(a) = (JL U M W]U V) W 
,dW 



WAu^W.u^W 



L 2 (fl) 



L 2 (T- 



A(u M W),u^W 



L 2 (Q) 



\ uM \ 2 g^ dx! + (K + n)(u^\u^W) L2i r E) 



({b + K + (j)«("'f, uMw ) L2 (r_) - Rc (wAu^W, u^W) 
Re(u {fM) AW,u M w) - 2Re (wVu^.u^VW) 

V J L 2 (Q.) V /L 2 (0) 



L 2 (fi) 



We calculate the last term in the right hand side of this equation by one more integration 
by parts. Using (2.9) and (3.2), we get 



-2Re(M / Vu (Al) ,w ( ' i) Viy) L2(n) = -- / VW 2 • V|u (/i) | 2 dx 

n 

= 1 f W 2JL\ u M\2 dxi + l f W 2 A\u^\ 2 dx 
Z J UX2 z J 



({b + K + LL)u (M) W,u (tl) W) 1 



r) + Re(u^W, WAu^) L2{n)) + ||WVu^||i a(n) . 
We substitute the obtained identity into the previous one and we find 



\\V(uMW)\\l 2m =J\uM\ 2 ^dx 1 + (K + ri( 



«00,u</0W) La( r.) 
/«Mll 2 - 



(3.38) 



-Re(uMAW,u^W) L2{n) + \\WVn- llL2{ny 
It follows from (3.31) and (3.35) that 

\(u^AW,u^W) L2m \^Ce^ 2 (K + ^ 2 \\ne(K + fi)\\\f\\l 2{ny 
Following the arguments of the proof of Lemma 3.5 in [1], one can check easily that 



\uM\^d Xl 

0x2 



^C(K + ri\\f\\l 2in) . 



Due to (3.36) wc have 

\(u^,u^W) L2irs) \^Ce(K + [i)\\ne(K + n)\\\f\\ L2m . 

We substitute the three last estimates and (3.31) into (3.38) and it completes the proof. 

□ 
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It is clear that the resolvent (H^ + i) 1 is analytic in /i and this is why 



\Wg> +1)- 1 - rf + i)~ 1 IU a (nHwJ(n) < <?/.. (3-39) 

This inequality, Lemma 3.6, the estimates (2.14), and (3.31) yield the estimate (2.12). 
The estimates (2.10) and (2.11) follow from (2.14), (3.31) and (3.39). 

4 Resolvent convergence: homogenized Dirichlet con- 
dition 

In this section we study the uniform resolvent convergence for the case of the homog- 
enized Dirichlet condition and prove Theorem 2.2. The case of periodic alternation 
■& £ (t) = t, a- = rj for the pure Laplacian T-L e = —A was proved in [4, Sec. 2]. Here we 
employ similar ideas but with changes required by the non-periodicity of the alternation 
and for the more general operator. 
Given / 6 Z,2(f2), let us denote 

u e := (H s - i) -1 /, "o := C^d - i) -1 /, v e ■= u e - u . 

The last function belongs to W-zotfl, T+ U 7e) and it is the generalized solution of the 
boundary value problem 



((iV + A) 2 + V-i)i 



in 0, 



v E = on T+ U 7e , [—faT ~ iA ' 2 + b ) Ve = a~ on Fe ' 

We multiply this equation by v e , integrate by parts over ft taking into consideration 
the boundary conditions, and get 



= v e (—+iA 2 )v e dx 1 + || (iV + A)i 



l£ 2 (fi) 



+ (Vv £ ,v £ ) L2(n) -i\\v s \ 



2 
L 2 (f2) 



|(iV + A)v £ \\ L2{n) + {Vv ei v e ) L2(n) - i||v e || L2(0 ) + (bv s ,v e ) L2{Te) - I 



fdu 

-= — ,f f 



9a; 2 



i 2 (r £ ) 
(4.1) 



Since u Q G W$ (ft) , the term ( ff 1 , v e ) is well-defined by the embedding of W| (11) 

into W^r.) and 



dup 

(j.r- 



■,Vi 



L2(T e ) 



<c 



du 



dx 2 



L 2 (T e 



MU 2 (r e ) < C||«o||w|(n)II^IU2(r e )- 



Here and till the end of the proof we denote by C inessential constants independent of 
e and /. The last estimate, (4.1), (3.37), and the identity v e = u e on T e yield 



IMIwjcn) ^ c (ll«e|li 2 (r e ) + IM|wf(n)IKIU 2 (r s ) 



(4.2) 



We employ again the rescaled variables £ defined in (3.9) and introduce one more 
auxiliary function 



Xv = x v(0 '■= Rem (ship + y sin p - sin rj) - £ 2 , P = 6 + 16, 

where the branches of the logarithm and the square root are fixed by the requirements 
in 1 = 0, yl = 1. This function was introduced in [22]. It was proven that it is harmonic 
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in the half-space £2 > 0, even and 7r-periodic in £1, decays exponentially as £2 ~> +00, 
and satisfies the boundary condition 



X n = In sin 77 

dx„ 

56 



on 7(77) := |J{£: |£i-7rj| < 77, 6 = 0}, 



-1 



on T( V )--OZi\j(r 1 ). 



In [4, Sec. 3] it was shown that the function X v is continuous in {£ : £2 ^ 0} and obeys 
the estimate 

\X V (0 K|ln sin r,| (4.3) 

uniformly for £2^0. 
Denote 

?7»(e) :— minjinf a+(e),inf aj(e)}. 



By the assumption (A3) we get 



V(e) < V*(e) < c 3 77(e) sC -. 



(4.4) 



In the same way as in Lemma 3.2 in [4], one can check easily that uA,( E )('£ 1 ) G 
VK 2 1 (ri,r + U 7 S ). We choose this product as the test function in the integral identity 
for u £ , 

i) £ [u £ ,u £ X Vt ] = (f,u e X n ,) La (a), 

or, by (2.3), 

((X7 + A)u E ,X v AiV + A)u £ ) L2m + ((iV + A)u £ ,iu £ VX Vt ) L2in) 

+ (Vu e ,X n ,U e )L 2 (n) + Q>U £ ,U £ ) L2 (r e ) -i(We)-X"f7.We)i 2 (fi) = (/, UeX v ,) L2 (Q), 

where X Vr = X I)t (-e _1 ). 

Employing the described properties of X v , we integrate by parts as follows, 

Re((iV + J 4)u e ,m e VX 7) ,) i2(n) = Rc(yu £l u £ VX nt ) L2{n] 

= - / VI,, • VK| 2 dx = -- \ue\ 2 -g^ dxi - / \u £ \ 2 AX Vt dx 



J- f tf' e \u £ \ 2 d Xl - J \u £ \ 2 AX Vr dx. 



2e 

r E n 

This identity, (2.5), and the real part of (4.5) imply 



he" 2 



j+2e(bu e ,X Vm u e ) L2 (r e ) =e / |u e | 2 AX,,, dx + 2e(/,X r; ,u £ ) L2(n) 

h 
- 2e((iV + i4)« ej Jf,.(iV + A)u £ ) L2{n) - 2e{Vu £ ,X, h u £ ) L2{n) . 

By the boundedness of b, V, and A it follows that 



Mli 2 (r e ) < Ce 



|uJ 2 AX„ da: 



+ |liiBin77,( e )|(||« B ||^i (n) + ||/||! 2( a)) • (4-6) 



Let us estimate the first term in the right hand side of the last inequality. 
In view of the harmonicity of X Vil and the definition (3.9) of £ we get 
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Hence, we can integrate by parts as follows, 
J \u E \ 2 AX n , dx = -2 J X^^-^u^ 2 dx - f X v ,-^-\u s \ 2 dx 



■&'" -&' 9 



2$'i 



<2 d& { V* 



7?" \ 2 \ d 2 \ 



and thus 



\u e \ 2 AX Vr dx 



sC C sup 



!#%. 



^ 2 oe 2 



sup 



dXr 



96 



- sup \X V , 

«2^0 



\\ u ^\\w^(n)- 

(4.7) 



Lemma 4.1. The estimate 



sup 

«2^0 



2 d X, h 



d& 



sup 
6>o 



dX 



'/. 



96 



^ C cos 77* 



holds true. 

Proof. Since the function X Vt is even and 7r-periodic in £1, we have 



sup 

?2^0 



.2 ^ -^»7» 



d& 



sup 

«2>0 



9X 



96 



= sup 

0^£i<t/2 
£2^0 



^ 2 a6 2 



sup 

«2^° 



dx. 



96 



(4.8) 



It follows from the definition of X Vt that 



dX r . 



96 



— Im 
-Im 



cosp 



ysin 2 p — sin 2 77 



- 1 = - Im 

cos 77* 



cosp 



/ 1 9" • 2" 

+ lysm p — sin 77 



V sin 2 p — sin 2 77 



ysin p — sin^ 77 



; n 2 ~ / cos p 
cos 77* 



fc" 1 ) 



As ^ 6 ^ 7r /2, 6^0) the fraction cosp ranges in the closure of the first quarter 
of the complex plane. In this quarter the function z 1— > z + \J z 2 — 1 has no zeroes. 
Moreover, as z tends to infinity within this quarter, the function z + y/z 2 — 1 tends to 
infinity, too. Hence, there exists a positive constant C such that 

\z + Vz 2 - 1| > C > uniformly in Re z ^ 0, lmz>0. 
This estimate implies 



cosp 



and therefore 



6 



cos 77* 



9X 



/ cos 2 p 
cos 2 77* 



^ C, 0^6^ tt/2, 6 ^ 0, 



96 



< 



C6 cos 77* 



|sin 2 p- sm 2 -!]^ 1 / 2 



, sC 6 < tt/2, 6 > 0, 



(4.9) 



where the constant C is independent of £ and 77* . By direct calculations we check that 

I sin 2 p — sin 2 77* I l ' 2 = I sin(p — 77*)sin(p + 77*)! 1 ' 2 

ma ma ( 4 -!0) 

= (sinh 2 6 +sin 2 (6 - ? 7*)r >inh 2 6 + sin 2 (6 +?7*)) , 

I sin 2 p- sin 2 77, I 1 / 2 ^ sinh6- (4-11) 
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Hence, by (4.9) 



dX r , 



56 



g C & cosr >* ^Ccosr]*, 0<£i<7r/2, 6>0. (4.12) 



sinh£ 2 

In the same way we estimate the first term in the right hand side of (4.8). We have 
d 2 X„ 



d£ 



cos 77* Re — 



sinp 



(sin p — sin 77* ) 3 / 2 



Thus, 



,a 2 x 



el-H^Wicos^Re — 



1 



^: 



(sin p — sin 77* )!/ 2 (sin p — sin 77*) 
1 



£2 cos 77* sin 77* Re — - 



^2 ^ j^g. 



03 



<- 



(sin p — sin 77*) ' 
£f cos 2 77* £f cos2 7 ?* sm 7 7* 



sin p — sin 77* | ^^Z 2 1 sin p — sin 77^ | | sin p — sin 77* | 3 / 2 
It is easy to see that 

I sin p + sin 77* | = sinh £2 + sin £1 + sin 77* + 2 sin £1 cosh £2 sin 77* Js sinh £2 
for < & < tt/2, £ 2 > 0. It also follows from (4.10) that 

I sin 2 p - sin 2 77, | 3/2 ^ ( sinh 2 £ 2 + sin 2 (& + 77)) 1/2 sinh 2 £ 2 

^ sinh £2 min (sin 77* , cos 77* ) Jj sinh £2 sin 77* cos 77* . 

We substitute two last estimates and (4.11) into (4.13), 

£2 cos 2 77* £! cos 77 



(4.13) 



2 d X Vr 



d& 



€ 



+ 



s^Ccost?*, 6^0, 0s^£i< 



sinh £2 sinh £2 
The obtained estimate, (4.12), and (4.8) complete the proof. 
The proven lemma and the estimates (4.3), (4.7) yield 



\u E \ 2 AX Vt , dx 



s^ C*(|msin77»(e)| +cos7/*(e))||w e || 2 V2l(0) . 



Together with (4.6) it implies 

IMl! 3 (T.) ^ C 4\ hls[lll h(t)\ + co sV*(£))(\K\\wi {n) + \\f\\l 2( n)) 
We observe that the functions uq and u e satisfy the estimates 

IKHw 2 2 (n) < C||/IU 2 (n), \\ u s\\wi(n) < C||/IU 2 (fi)- 
Due to (4.4) we have 

I In sin 77, (e) I ^ C||/|| L2 ( ), cos77*(e) ^ cos?7(e). 
It follows from last five estimates and (4.2) that 

IKII^cn) < ^ 1/2 (l Insin 77(e)! + cos 77(e)) 1/2 1| /||i 2(n) 
and it proves Theorem 2.2. 



□ 
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5 Analysis of the operator H e (r) 

In this section we prove Theorem 2.4. 

We state now two auxiliary lemmas proved in [4] and [1] that we give below for the 
reader's convenience. 

Lemma 5.1. Let |r| < 1 — *c, where x € (0, 1), and for a given function f G L 2 (£l e ) 
we let 

u e = (n e (r) - ^) /, / g L 2 (n E ). 

Assume f G £^r . Then 

£ £ 

\\U e \\L 2 (n e ) < — jToll/II^Cfie)) ll V ^ell-L 2 (n e ) ^ 7T~ ll/IU 2 (n e ) • 

For any u G W 2 per {£l e ,T +) and |r| ^ 1 — x, we have the following inequalities 



. d t 
i^ I 

ax\ e 



L 2 (n e ) 



^l|w|L 2 (n e) >* 



On 



dx\ 



L 2 (fi e ) 



da 



dx 2 



L 2 (O e ) " 2 11 ' 



ia(«)- 



Lemma 5.2. J/F G L 2 (0,7r), ifcen Kfi^-FXO)! < 5\\F\\ L2{0 ^. 
Let / G L 2 (n e ) and / = F £ + fr, where F £ G £ e , ft G *£, 

"2" 

1 



F e {x 2 ) 



TVS 



1/2 



/ £ (x)dxi, fe(x):=f(x)-F e (x 2 ), 



(5.1) 



^II^IIW) + 11/, lli 2(ne) = 11/111,(0,), U e := ( H s {t) - ? 



F K . 



It follows from Lemma 5.1 that 



^e(r)--. /-J7 8 



L 2 (n £ ) 



^W -72 /e X 



(5.2) 



L 2 (fi e ) 



^- 1/2 ||/IU 2( n s) . 

(5.3) 



Denote 



C/i^ := Q^F e , V e {x) := C7 e (x) - U^\x) - U^(0)W(x,e,n) X (x 2 ). (5.4) 

We remind that the function \ was introduced after the equation (3.21). The boundary 
conditions (3.2) and the definition of U £ imply that V e G W 2per (Q, e ,Y + U %). 

It follows from the definition of U e and Ue that they satisfy the integral equations 

^2 



Ur){U £ ,V £ ] - -(U e ,V s ) L2(ne) = (F e ,V e ) La(flt) , 



(5.5) 



% [U £ ,V e ]- ^{U s ,V s {x ir )) L2{0n) = {F e ,V s (x u -)) L ^ y 
We integrate the last equation over x\ G (— £7r/2,£7r/2) and use an obvious relation 



ctei e/ £ ' \ dx i e 



L 2 (n e ) 



9^1 £ 



i 2 (n e ) 



= ^(^),T4) L2(0e) . 
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It leads us to the identity 



Vw (1A-1W ^ ' 



dx\ el V dxi e I I T re, \ \ 8x2 ' dx2 I 



- ^([/^),y £ ) L2(0e) + (6 + /i + M )(C/^,K) L2( f e) - (^,K) i2 (a £ )- 



^ 2 
We take the difference between the last identity and (5.5), 



ace 

U e -U^ = V £ + U^(0)W X , 



we get 

r 2 



Ur)[V e ,V £ }- ^\\V e \\l 2{ne) =(K + ^)(U^,V e ) L2{te) 

- uM(0)t>e(T)[Wx,V e ] - ^(W X ,V e ) La{ne) . 
Integrating by parts and employing (3.2), we obtain 

Ur)[W X ,V e }-^(W X ,V £ ) L2m 

it ( dWy \ \t ( dV F \ 

=(v^,w,)„ , ) - _ {^v,) lmi + 7 (^) MCU + W V. W ., 

— 9Wy 2ir / <9K\ 

G ' x 2 t \ G ' x l/L 2 (r! e ) 



= - (AW X , K)i 2 (n e ) + — ( W X , ^-) + (bW + K + n, V s ) 

e \ ax i/L-,(n r ) 

The last two equations imply 



L 2 (r e y 
L 2 (n e ) 



be(T)\y e ,V e ] ' ^\\Ve\\l 2( n e ) = -U^(0)(bW, V,)^ f) 

+ U^(0)(AW X ,V s ) L2{ne) - — J7J">(0) (w X ,^ 

e \ axiy i2 (n e ) 



(5.6) 



Lemma 5.2 and the inequalities (5.2) give the estimate for U £ (0), 

Ic/^WKs^e)- 1 / 2 !!/!!^^). (5.7) 

In the same way as in Theorem 2.3 in [1], we get the estimates for Wx and AWx, 

\\W X \\L 2 (n e ) < \\W\\ L2( a s ) ^ Ce 2 (K + M ), 
IIAWxIU^o^^iir + ^^ + e-- 1 ). 

Using the last three inequalities, we estimate the second and the third term in the right 
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hand side of (5.6) as 



uM(0)(AW X ,V e ) L2( n e) - ^fu^(0) (w X ,j^ 

< 5(^r i/2 n/iu 2( o) mia^xiIl 2 (o e )||kiil 2 (o e ) 



L 2 (n c ) 



+ 2 £ - 1 ||^ X || L2(ne) 



8V R 



dxi 



(5.9) 



l 2 (q c ) 



<> 



dV e 



dxi 



L 2 (n e ) 



+ T^\\^\\Un B} + C^e(K + ^\\f\\i 2(ne) 
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It remains to estimate the first term in the right hand side of (5.6). We employ the 
obvious inequality 

||V e |lr_rf.^C 



' £ "L2(r e 



8x2 



L 2 (Q e 



(3.25), and (5.7) to obtain 



\uM(0)(bW,V s ) L ,, J ^ C\U^(0)\\\bW\\ L , t jV 6 



>L a (T.)l ^l"e V»;ill""MlL 2 (r e )ll^llL 2 (r e ) 

8V. 



^Ce-^Wfh^) 



3X2 



1 
< - 

2 



8V F 



dx 2 



L 2 (Qe) 

2 



w K<t.^'< ) + \\ w K<t.rt*'< ) 



L 2 (n e ) 



+ Ce(K + „y\lne(K + ^\\f\\i 2iney 



We substitute this estimate and (5.9) into (5.6), and by (5.1) it follows 



dV, 



dx\ 



L 2 (a e ) 



1 

+ 2 



8V F 



8X2 



-UK 



L 2 (fi e ) 



L 2 {n e ) 



€ 



8V F 



8x\ 



L 2 (n E 



+ Till^llLtoj + Ce{K + tf (x" 1 + | \ne{K + M )| 2 ) ||/||| 2(0e) . 
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Therefore, 

\\V 6 \\ L2{ n E ) ^ Ce l ' 2 {K + (i)^ 2 + | lns(K + M)|)||/IU a( n.)- 

Together with the definition (5.4) of V 6 , (5.7), and (5.8) it yields 

\\U S - U^ || L2(i2e) = \\V e + U^(0)W X \\L 2 (n E ) 

< Ce^iK + Af )(^ 1/2 + I lne(K + M )| + e)||/|| L2(f , e) . 

We combine this inequality with (5.3) and it completes the proof of (2.19). 

The proof of the asymptotics (2.20), (2.21) is similar to that of Theorem 2.4 in [1]; 
it is enough to use the proven asymptotics (2.19) instead of Theorem 2.3 in [1]. 

6 Bottom of the essential spectrum 

In this section we prove Theorem 2.5. 

We begin by proving the identity (2.4). We first observe that the form \) f (t) associ- 
ated with the operator T-L s (r) can be estimated from below as 



t) e (r)[u,u] ^ 



. d t 

1 n I 

oxi e 



L 2 (n e ) 



On 



8xn 



L 2 (fi e ) 



i£ 2 (r_)> 
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where the form in the right hand side is on VPg (0 E ,f + ). Hence, by the minimax 
principle we have 



Ai(e,r) ^ inf 

H"Hl. 2 (fi)=l 



d t 

1 n I 

oxi e 



L 2 (n e ) 



On 



dx? 



L 2 (O e ) 



lia(r-) 



(6.1) 

The infimum in the right hand side of this estimate is the lowest eigenvalue of the 
operator 



~h in L2{ne 



d_ 
dx\ 



with periodic boundary condition on the lateral boundaries of f2 e , the Dirichlct condition 
on r + , and the Robin condition 



b\u = on I\_. 



_d_ 

dx 2 

We find this eigenvalue by using the separation of variables and substitute it into (6.1), 



Ai(e,T)>-j +T\ 

where T is the smallest nonnegative root of the equation 

Tcos7rr + &sin7rT = 0. 
It is clear that the first root of (2.21) can be estimated uniformly in small /i as 



(6.2) 



^ \J k\{n) < c, c = const. 
Thus, by (6.2) for |r| ^ (c+l)e, r £ [-1, 1) we have 



Since by (2.20) 



Ai(e,r) > (VAi^ + l) 2 > Ai(/*) + 1. 



X 1 (s,0) = A 1 ( l x)+O(e 1 / 2 ), e- 
we conclude that for sufficiently small e we have 

inf Ai(r, e) = inf Ai(r, e). 

re[-l,l) |tK(c+1)£ 

Relation (2.4) is now proven by the arguments used in [4, Sec. 5] and based on Temple 
inequalities. 

The rest of the section is devoted to construction of the asymptotic expansion for 
Ai(0,e). Here we employ the approach suggested in [14], [21], [23], [24]. 

We write the boundary value problem for the eigenvalue Ai(0, e) and the associated 
cigenfunction tp(x,e), 



ijj = on T+ U 7, 



AV = Ai(0,e)^ iu ih, 
d 



b\ip = Q on r £ , 



(6.3) 



and on the lateral boundaries of Q e the periodic boundary conditions are assumed. The 
asymptotics for Ai(0,e) is constructed as 

\ 1 (0,e)=A(jt,e), (6.4) 

where the function A(/x, e) is to be determined. By (2.20) it should satisfy the identity 

A(e,n)=A 1 {n) + 0(e 1 > 2 {K + n)), e -> +0, y, ->■ +0. (6.5) 
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The asymptotics for ip is constructed by the combination of the boundary layer 
method [39] and the method of matching asymptotic expansions [38] . It is sought as 
a sum of an outer expansion, a boundary layer, and the inner expansion. The outer 
expansion is defined as 



V>f(x, A) = sin y/A(ji,e)(x2 - tt). (6.6) 

This function is periodic w.r.t. x\ and satisfies the equation and the boundary condition 
on r + in (6.3) no matter how the function A looks like. 

The boundary layer is constructed in terms of the rescaled variables £, and we 
denote it as ip^ l (^ £*)• The main idea of using the boundary layer is to satisfy the 
required boundary condition on T e and thus 

JL-bj(J? + i?) = on f £ . 

We substitute (6.6) into this condition, rewrite it in the variables £, and pass to the 
limit as r\ — > +0. It leads us to one more boundary condition, 

^--e6J^ 1 = -e(\/Acos\/A7r + &sin\/A7r) on f°, (6.7) 

f°:={e:0<|e 1 |<|,6>0}. 

To obtain the equation for rpf-, we substitute it and (6.4) into (6.3) and pass to the 
variables £. It yields the boundary value problem 

-A ? ^ = e 2 A^, cen, n:={M6l<|, 6>o}, (6.8) 

and on the lateral boundaries of II the periodic boundary condition is imposed. 

By 03 we denote the space of 7r-periodic, even in £1 functions belonging to C°° (II\{0}) 
and exponentially decaying as £2 — > +00, together with all their derivatives, uniformly 
in £1. It is easy to check that the function X introduced in (3.12) belongs to 03. We 
state the next lemma that was proven in [1]. 

Lemma 6.1. The function X can be represented as the series 

00 1 
X(£) = -^-e- 2 "« 2 cos2<i, (6.9) 

n=i n 

which converges in 1*2(11) and in C fe (II n {£ : £ ^ R}) for each k ^ 0, R > 0. 
Lemma 6.2. For small fj the problem 

-A/:Z- (3 2 Z = in II, ( — - eb\ Z = on f°, (6.10) 

has an even periodic in £1 solution in 03. This solution and all its derivatives w.r.t. £ 
decay exponentially as £ 2 — ► +°o uniformly in £i, e, and (3. The differ entiable asymp- 
totics 

Z{£, e, (3) =c^ (in |£| + In 2 + 9(eb, (3 2 ) - (1 + ebfe + ehfr [tp - |)) 

+ 0(|£| 2 ln|£|), C^O, 

holds true uniformly in fJ and e, where ip is the polar angle associated with £, and the 
function 6 is defined in (2.22). The function Z is bounded in £2(1!) uniformly in e and 
j3. The function 6{t\,t2) is jointly holomorphic in t\ and t^. 
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Proof. We make the change 

/ +00 



Z(£,e,/3)=e^ 



in the problem (6.10), 



Z(Z,e,P)+X(S)-ebS2X-eb f X(£ u t)dt\ (6.12) 



-A-2eb^-£ 2 b 2 -f3 2 )Z = F in n, %— = on L°, (6.13) 

(+00 \ 

&X+ J Xfo,t)dt , 
& / 

and on the lateral boundaries of II we have the periodic boundary conditions. 

Let 03 be the orthogonal complement in £2(11) to the set of the functions <fi = 0(6) 
belonging to £2(11). By 2U we denote the Hilbert space of functions in Wf (II) n 03 
satisfying periodic boundary conditions on the lateral boundaries of II, the Neumann 
boundary condition on L°. In 03 we introduce the operator B in acting as —A and on 
the domain 233. This operator is self-adjoint. 

It was shown in the proof of Lemma 5.2 in [1] that B ^ and hence the inverse 
exists and is bounded. It follows that 2? _1 is also bounded as an operator from 03 into 
W2 (II) n 03. Hence, for sufficiently small e and /3, we have 

B . 2eS |_ _ sV . p) - . b-. (, . ( 2eS |_ + eV + ^) B - 

where we have beared in mind that the operator ■§- maps W^ 1 (II) n 03 into 03. 
Thus, the solution to the problem (6.13) reads as 

a \ _1 

B - 2eb- e 2 b 2 - /3 2 ) F. 

96 / 

By the standard smoothness improving theorems we conclude that Z <G C°°(II \ {0}). 
The function Z can be found from (6.12), 



B - 2eb^- - e 2 b 2 - ($ 2 



-1 + °° \ 

F + X-eb&X-eb f Xfa,t)dt\ 

i I 



We can obtain one more representation for Z by the separation of variables. In order 
to do it, we construct Z as 

Z(£,e,l3) = Z(Z,e,p) + e^X(Z) 
and then we separate the variables for Z . It implies 

Z(t,e,t3) = -J2 lA l m -cos2n6. (6.14) 

~ sj\v? -(3 2 +eb 

As in the proof of Lemma 5.1 in [1], we check that this series converges in £2(11) and in 
C fc (TTn {£ : 6 ^ R}) f° r eacn k ^ 0, R > 0. It implies that the function Z and all its 
derivatives decay exponentially as 6 ~~ ^ +°° uniformly in £1, e, and /3. Hence, Z <G 03 
and for sufficiently small e and /3 
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~ ^ (V4n 2 - 2 + eb) V4n 2 - /3 2 ^ ^ (4n 2 - 1)3/ 2 ' 

Hence, the function Z is bounded in L 2 (Ji) uniformly in e and 0. 
By analogy with Lemma 3.2 in [42] one can prove that 

Z(£, £ ,/3) = Z(O,£,/3) + 0(|£| 2 ln|£|), £^ 0, (6.15) 

uniformly in small e and 0. Let us calculate Z(O,e,0). 

We employ the asymptotics (3.13) with j = and (6.15) and we integrate by parts 
as follows, 

JxAZd^lun J {z^-X^\ds = ,Z%e,0). 

n iei=*> ^ ' 

«2>0 

Hence, by (6.12) and the equation in (6.13) we get 

Z(0, e,0)=-±-Jx (( 2£b -^T + **& + P 2 )z + f) dti 
n 

-^/ JC (( 2£6 i +£V+ ^) e "* z - &6 I)*- 

n 
We substitute the series (6.9), (6.14) into the obtained equation, 

Z(0 e 0) = Vl /" f 2£& V^^ 1 + ^ ~ ^ c ^ ( ,/^^ + 2n + e b ) 

^i n J \ V^ 2 -0 2 +eb 

n—i 

~2ebc~ 4n ^) d£ 2 



_^yI + M y : (6.i6) 

2 „ti« 2 J^l n ( ^4n 2 - /? 2 + 2n + eft) 

-(/3 2 +e v)$: 



ri n( y/An 2 - 2 + eb)(y/4n 2 - 2 + 2n + eb) 

= -^ + e( sb ,n 

In view of (6.12), to prove (6.11), it is required to study also the behavior of the 

function J X(£±,t) dt as £ — > 0. We have 

€2 

+ OO +OO £2 

j Xfa,t)dt= J X(^t)dt- [x(£ u t)dt. (6.17) 

£2 

It follows directly from the definition (3.12) of X and (3.11) that 

+00 +00 



^ j x(ti,t)<tt = - J ^x(£ 1 ,t)dt = -i, a 7^0, 



+00 +00 

/'(40*-s- /f(4<) 
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and thus 

+00 







The asymptotics (3.13) implies 



/"x(a,t)dt-61n|el + (ln2-lK 2 +aarctan^ + 0(|e| 2 ), £ -► 0. 


The two last equations, (6.17), and (3.13) yield 

+00 

6^(0 + / *(&, *) dt = ~ + 6 (| - *>) - 6 + o(\e), t -> 0. 

This identity, (6.12), and (6.16) lead us to (6.11). The function 9 is jointly holomorphic 
in t\ and £2 by the first Weierstrass theorem since the terms of series (2.22) are jointly 
holomorphic in t\ and £2 and these series converges uniformly in small t\ and £2- D 

The proven lemma allows us to construct the needed solution to the problem (6.8), 
(6.7). Namely, we have 

$?'(£, e, A) = e(%/Acos VAtt + 6 sin y/Lir)Z(£, e, sVX). 

By the definition (6.6) of ip° x , it has the asymptotics 

■0° x (x, A) = - sin VXtt + VXx 2 cos VAtt + 0{x\) 

= e bX2 (-sinVA7r+ (\/Acos\/A7r + 6 sin \/A7r) x 2 ) , x 2 -> +0. 

The asymptotics of ip^ 1 as £ — > can be easily found by (6.11), and, as £ — >■ 0, we have 

^ X (*, A) + $1 ( £ E> A) =e e^ 2 ^in (?) £j A) + £r? £in (?) ^ A) ) 

+ 0( £ 7 ? 2 k«| 2 (|lnk||+ £ - 1 (X + / i))), 

*[, n (<?, e, A) =e(a/X cos VAtt + b sin VAtt) (In |?| + In 2) 

+ (V~A cos VItt + b sin VAtt) ( e 0(e&, £ 2 A) - (K + /i)" 1 ) (6.18) 

— sin vAtt, 
*i n (<T,e,A) =e6(VAcos\/A7r + 6sinVA7r) ( ?1 (</> - |) -ft) , (6-19) 

where <j = £?7 _1 . In accordance with the method of matching asymptotic expansions we 
construct the inner expansion for ip in a vicinity of the point x = as 

C'^A) = e e ^(C( ? , £ ,A) + £ r 7 ^( ? , £ ,A)), (6.20) 

where the functions ipi, i = 0, 1, must behave at infinity as 

C^,e, A) = *j n (,, £ , A) + 0(kP), ?^oo. (6.21) 

We substitute the ansatz (6.20) into the boundary value problem (6.2) and equate the 
coefficients at the like powers of rj. It implies the boundary value problems for ip™, 

A^ n = 0, ? 2 >0, C = 0, CS7 1 , ^ = 0, ? er\ (6.22) 

A^ = -2b^, , 2 >0, ^ n = 0, C6f, if^-O, seT 1 , (6.23) 

C?2 C^2 
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where 7 1 := {? : |ft| < 1, ft = 0}, f 1 := Oft \ 7 1 . The problem (6.22) has the only 
bounded at infinity solution which is trivial. The solution behaving at infinity as In |c| 
is also unique and it is the function Y introduced in (3.18). Hence, 

$> n (ft£,A) =£(VAcosV / A7r + 6sinV / A7r)F(?). 

The asymptotics of this function as cr — > 00 can be found by using (3.19). Comparing 
this asymptotics with (6.18), (6.21), we arrive at the equation 

(\/Acos\/A7r + 6sin\/A7r) (e6»(£&,£ 2 A) - (K + /j.y 1 ) -shiv^ = 0, 

which can be rewritten as (2.25). In the same way how similar equation (2.15) was 
studied in [1], one can prove easily that the equation (2.25) has the unique root satis- 
fying (6.5). This root is jointly holomorphic in e and /1 and it can be represented as 
the convergent series (2.27). To calculate explicitly its coefficients, it is sufficient to 
substitute this series into the equation (2.25) and expand it into the Taylor series w.r.t. 
£, set the coefficients at the powers of e equal to zero, and solve the equations obtained. 
Exactly in this way, one can check the formulas (2.28). 

To solve the problem (6.23), (6.19), we need an auxiliary lemma. 

Lemma 6.3. The problem 

BY BY 

AY! = -2— , ft>0, r 1= 0, cef, Tr 2 ^ ' ?efl , (6.24) 

Oft Oft 

has a solution with the differentiable asymptotics 

F 1 (?)=ft(^-|)-?2 + c+0(M- 1 ), S-+00, (6.25) 

where c is a constant. This solution belongs to W$ \{Q) H C°°({<; : ft ^ 0} \ {<; : ft = 
±1, ft = 0}), where Q is any bounded subdomain of {s : ft ^ 0}. 



Proof. Denote 



Fi (1) (?) := ft Imln (z + \J ' z 2 - l) - Im \j z 2 -1 - -ft. 



It is easy to check that this function solves the equation in (6.24), satisfies the boundary 
condition on f 1 in (6.24) and the asymptotics (6.25), and belongs to W\(ff) fl C°°({<; : 
ft > 0} \ {^ : ft = ±1, ft = 0}). 

Let xi = Xi(t) be an infinitely differentiable cut-off function with values in [0, 1], 
being one as t > 3 and vanishing as t < 2. We construct the solution to the problem 
(6.24), (6.25) as 

r 1 (0 = r 1 (1) (0 + xi(kl)r 1 (2) W- 

The function Y± should solve the problem 

Ay/ 2 ' = -2(1 - X i)|^ - 2V X i ■ vyW - y^Axi, ft > o, 

Oft 
Oft 



n (2) (0 = o(kl), 



00. 



The right hand side in the equation of this problem is compactly supported due to the 
definition of xi and belongs to X 2 ({c : ft > 0})nC*°°({<; : ft ^ 0}\{? : ft = ±1, ft = 0}). 
At the next step we make the Kelvin transform and pass to new variables 

/ s ft ft + 1 

2/ = (2/1,2/2), 2/1 = ^—7 — -TT2' 2/2 = 



ft 2 + (ft + i) 2 ' " ft 2 + (ft + i) 2 

30 



(2) 

Under this transform the boundary value problem for Y± casts into the boundary 
value problem for the Poisson equation in the disk {y : y\ + (j/2 — 1/2) 2 < 1/2} with 
a combination of the homogeneous Dirichlet and Neumann condition. The right hand 
side in this equation is compactly supported with the support separated from zero and 
it belongs to Li. Hence, we can apply standard theory of generalized solutions to 
elliptic boundary value problems, see, for instance, [41, Ch. IV]. In accordance with this 
theory, the obtained boundary value problem is solvable in W^dy : y\ + (2/2 — 1/2) 2 < 
1/2}). By the smoothness improving theorems, the solution to this problem is infinitely 
diffcrcntiable in a vicinity of zero. Returning back to the problem (6.26), we conclude 
that it is solvable in W^Q) within the class of bounded at infinity functions. Hence, 
the problem (6.24) is solvable and has the solution with the asymptotics (6.25). By the 
smoothness improving theorems the solution belongs W^Q) H C°°{{<, : <?2 ^ 0} \ {<; : 
ft=±l, <^2 = 0}). □ 

In view of the proven lemma the solution to the problem (6.23), (6.19) is 

V4 n (?,e,A) = £&(\/Acos\/A7r + 6sin\/A7r)Yi(<;). 

The formal construction of the asymptotics is completed. 
Denote 

* e (x) := (4 cx (x, A) + x(*2)4 bl (^ A)) (1 - X ) + x{\^ 1/2 W n {^ A), 

where A = A(e, 11(e)) is the root to the equation (2.25). The proof of the next lemma 
is analogous to that of Lemma 5.3 in [1] and is based on direct calculations. 

Lemma 6.4. The function ty E belongs to C°°(f2 e \ {x : x\ = ±£77, xi = 0}) and to the 
domain of % e (0) . It satisfies the convergence 



|| * e - sin VMO) (a- 2 - 7r)|| i2(ne) ->• 0, s -> +0, 

and solves the equation 

(7i e (0)-A)* e = ft ej 

where for the function h e G £2(^5) the uniform in e estimate 

HMi a (n.) < C({K + /J )c- 2 ^ 1 + ^{K + /*) + e^W'^K + m) 1 / 2 ) 

holds true. 

Employing this lemma and proceeding as in [1, Sec. 5], one can easily check (2.24). 
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